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The parton content of the nucleon is explored within a meson-cloud model developed 
to derive light-cone wave functions for the physical nucleon. The model is here applied 
to study electromagnetic form factors, distribution amplitudes and nucleon-to-meson 
transition distribution amplitudes. 
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1. Introduction 

Schematically, in light-front quantization a nucleon state is described by the follow- 
ing superposition 

\n) = *P( 3 q) \qqq) + i>{z q +i g )\qqqg) + ^(M+qq) \qqqqq) + ■■■, (i) 

where in the light-cone gauge A + = the valence three-quark light-cone wave 
function (LCWF) ^(Sq) involves six independent amplitudes corresponding to dif- 
ferent combinations of quark orbital angular momentum and helicity, and the Fock 
component ^(39+13) with three quarks plus one gluon involves 126 independent 
amplitudes.^ Adding a pair of sea quarks into the valence component to build the 
amplitude ip(3 q + q q) leads to an even more complicated LCWF. 

To probe the parton content of the nucleon suitable models have to be invented to 
give explicit expressions for the LCWFs. Here, LCWFs are derived in a meson-cloud 
modeP^ 1 where the physical nucleon is pictured to be part of the time a bare nucleon 
and part of the time a baryon-meson system as a consequence of the spontaneously- 
broken chiral symmetry. Configurations are considered where the baryon can be a 
nucleon or a A and the meson can be a pion as well as a vector meson such as the p or 
the lo. The meson cloud will manifest itself as an extension of the charge distribution 
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of protons and neutrons, which should be observable in the electromagnetic form 
factors at relatively small values of Q 2 . The valence quark content of the LCWF will 
determine the distribution amplitudes (DAs) driving the soft contribution to the 
high Q 2 behavior of form factors. Further information on the LCWF can be obtained 
by studying transition distribution amplitudes (TDAs) occurring in the description, 
e.g., of hard exclusive backward pion electroproductior3 or in associated production 
of a pion and a high-Q 2 dilepton pair in pp annihilation.^ 

2. Light-cone wave functions and the meson-cloud model of the 
nucleon 

Following Refs. [2j [3] and references therein, in the one-meson approximation of 
the meson-cloud model the quantum state of the physical nucleon (N), with four- 
momentum p^ N = (p~^,p~^, Pjv_l) = (pJ[,pn) and helicity A, can be written as 

\p N ,X;N) =VZ\p N ,X;N)+j: f ^L—L=T1 ^ BM \VM 

x \yp+ , k ± + yp N± ,X'; B) |(1 - y)p+, -k ± + (1 - y)p N± ,X"; M), (2) 

where the function (f>y^y,' BM \y is the probability amplitude to find a physical 
nucleon with helicity A in a state consisting of a virtual baryon B and a virtual me- 
son M, with the baryon (meson) having helicity A' (A"), longitudinal momentum 
fraction y (1 — y) and transverse momentum k^ (— kj_). From the normalization 
condition of the nucleon state one obtains the following condition on the renormal- 
ization factor Z 

1 = Z+Y,Pbm/n, (3) 

B,M 

with 

f dyd 2 kj_ I ,l/2(N,BM) i , m2 fA \ 

P B M/N = J 2(27r)3 2-, I^A'A" ■ ( 4 ) 

Adopting a light-cone constituent quark model for the hadron states \ph,X;H) of 
the bare nucleon (N = 3) and baryon-meson (N = 5) components, and using the 
minimal Fock-state wave function in the light-cone formalism, one has 

-s/flSfrM 1 -!-) Sw) *"*KfH 

N 

X ({ x ii kiXj Xi,Ti}i—l t ,,, t jsr) J^J \XiPjj, Pi±,Xi,Ti), (5) 

i=l 

H \ f] 

where \E f A ' Ul ({xi, ki±; Aj, Ti}^!,..^) is the momentum LCWF which gives the prob- 
ability amplitude for finding in the hadron N partons (N — 3 and N — 2 for the 
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baryon and meson case, respectively) with momenta (xipjj,Pi± = kjj_ + XiPn±), 
and spin and isospin variables Aj and Tj, respectively. 

H f f] 

The wave function ^ A can be obtained by transforming the ordinary equal- 
time (instant-form) wave function \1>^'' C ' in the rest frame into that in the light-front 
dynamics, by taking into account relativistic effects such as the Melosh-Wigner 
rotation. The instant-form wave function is here constructed as the product of 
a momentum wave function, which is spherically symmetric and invariant under 
permutations, and a spin-isospin wave function, which is uniquely determined by 
SU(6)-symmetry requirements. 

In the case of the nucleon, the momentum wavefunction of Ref. |6| has been 
adopted. It reads 

N' 

- raw (6) 

with N' a normalization factor and Mq = Ylj—i w i — y/ 771,2 + k 2 ^ ne mass of 

the noninteracting three-quark system. In Eq. ©, the scale (3, the parameter 7 for 
the power-law behaviour, and the quark mass m are taken as free parameters. They 
have been fixecP to fit the magnetic moments pP and p n , the proton axial coupling 
constant qa = Ga(0), the Sachs form factors G 7 ^ at Q 2 = 0.15 GeV 2 , and G P E 
and G P M at Q 2 — 0.15 and 0.45 GeV 2 . A 5% uncertainty was allowed in the fitting 
procedure. A rather good fit is obtained in the proton case in the whole range of 
available data both for the Sachs form factors and the axial form factor Ga, while 
in the neutron case the fit is less satisfactory. In any case, the contribution from the 
meson cloud is smooth and mainly significant for Q 2 < 0.5 GeV 2 with a maximum 
at Q 2 = 0. The fitted values are 7 = 3.21, (3 = 0.489 GeV and to = 0.264 GeV. 
These values differ from the original set of parameters (7 = 3.5, (3 = 0.607 GeV and 
to = 0.263 GeVjSl, defined in order to reproduce fi p and /i" only. 

The canonical wavefunction of the pion is taken from Ref. |7| and reads 

^• [c] (ki,k 2 ) - ^/l aV2 exp[-fc 2 /(2a 2 )], (7) 

with k = ki = — k2, and the two parameters a = 0.3659 GeV and m q = 0.22 GeV 
fitted to the pion form factor data. The wavefunction of the p differs from the pion 
only in the spin component, with the rest-frame spin states of the qq pair coupled 
to J = 1 instead of J = 0. Similarly, the uj is described by the same spin and 
momentum wavefunction as the p, but with the isospin component corresponding 
to a singlet state. 

With such ingredients, the probabilities of the bare nucleon state, the (pir) and 
(An) configurations are 90%, 5.1% and 3.4%, respectively. 

3. Nucleon charge and magnetization densities 

It is known that the presence of a small admixture (1-2%) of mixed-symmetry 
S'-wave components greatly improves the result for the neutron electric form fac- 
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tor. Following Ref. [8] the mixed-symmetry S'-wave component is assumed to be 
represented by an appropriate combination of mixed-symmetry spin-isospin wave- 
functions with two radial wavefunctions of mixed symmetry. The consequences of 
including a small percentage of such mixed-symmetric contribution to the neutron 
electric form factor are illustrated in Fig. [TJ Even a percentage as small as 1% is 
able to produce a quite good result compared to data. In particular, including also a 
mixed-symmetry S'-wave component the neutron charge radius becomes quite close 
to the experimental value (r 2 = —0.112 fm 2 ). As anticipated in Ref. [8] the same 
calculation leaves the other nucleon form factors almost unaffected. 




Fig. 1. The electric form factor of the neutron. Solid line for the bare neutron wavefunction with 
SU(6) symmetry, dashed (dotted) line including 1% (2%) mixed-symmetry 5-statc. 




1 2 3 1 2 3 

b (fm) b (fin) 



Fig. 2. The proton and neutron charge densities as a function of the impact parameter b. Solid 
lines for a permutationally symmetric momentum wave function, dashed lines including mixed- 
symmetry components. 



October 7, 2009 8:34 



proc'trento 



Probing the parton content of the nucleon 5 

The radial distributions of the nucleon charge (fi c h) an d magnetization (p rn ) 
have been estimated for a long time by taking the Fourier-Bessel transform of the 
Sachs form factors in the Breit frame and neglecting relativistic corrections. The 
problem of unambiguously determining the charge density can be solved by looking 
at the charge density p(b) of partons in the transverse (impact parameter) plane 
as a function of the distance b relative to the transverse center of longitudinal 
momentum!^ This is possible because in the transverse plane boosts are purely 
kinematical, i.e. in the light-front framework they form a Galilei subgroup of the 
Poincare group. Then p(b) is the two-dimensional Fourier transform of the Dirac 
form factor F\: 

1 f°° 

P(b) = 7r dQQMQ^F.iQ 2 ), (8) 
27T Jo 

where Jo is a cylindrical Bessel function. 

The corresponding charge densities for the proton and the neutron are plotted 
in Fig. [21 They are concentrated at low values of b with a positive peak for the 
proton and a negative peak for the neutron, in agreement with the phenomenological 
analysis of Ref. 03 

These nucleon charge densities can be related to quark transverse distributions. 
Assuming that only up and down quarks are in the nucleon and invoking isospin 
symmetry, u(b) and d(b) can be obtained using 

u(b) = (f{b) + \p n {b), d(b) = pP(b) + 2p n (b). (9) 




b I fin I b lfm> 

Fig. 3. Transverse distribution of up and down quarks in a longitudinally polarized proton as a 
function of the impact parameter b. Left panels: the axial contributions Au and Ad for up and 
down quarks, respectively. Right panels: total contribution for quarks polarized in the longitudinal 
direction, either parallel (solid lines) or antiparallel (dashed lines) to the proton helicity. 
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Making additional use of the axial form factor, one can derive the probability 
p q (b, A, A) to find a quark with transverse position b and light-cone helicity A (= ±1) 
in the nucleon with longitudinal polarization A (= ±l)j^i.e. 

pP(J>,\&) = ±\pP(b)+\AAq(b)], (10) 

where p q {b) is defined in Eq. © and Aq(b) is the Fourier transform of G q A (Q 2 ). 
Assuming a positive proton helicity (A = 1) the resulting probability is shown in 
the right panels of Fig.[3l The axial contributions Au(b) and Ad(b) for up and down 
quarks (left panels), respectively, have opposite sign. When suitably combined with 
the corresponding transverse distributions u(b) and d(b) the positive helicity up 
quarks in the proton turn out to be preferentially aligned with the proton helicity, 
while the opposite occurs for down quarks. This result corresponds to that shown in 
Fig. 7 of Ref. [TUI where quite a different radial distribution of the axially symmetric 
spin density was presented for up and down quarks in the transverse plane. 



4. Distribution amplitudes 

The proton DAs are derived from the following proton-to-vacuum matrix elements 
of trilocal operators built of quark and gluon fieldJ^ 

D a,p n = F ^(0|e ufc u^(zin)[zi;z o ]i'iM^ (z 2 n)[z 2 ; z ]j>jd* (z 3 n)[z 3 ; z ]k'k\p(p P , A))) , 

(11) 

where T means Fourier transform, \p(p p , A)) denotes the proton state with momen- 
tum p p (p 2 — M 2 ) and helicity A; u, d are the field operators for up and down 
quarks, respectively; the Greek letters a, (3 and 7 stand for Dirac indices, while the 
Latin letters i, j and k refer to color; n is an arbitrary light-like vector (n 2 = 0) 
and Zi are real numbers that specify quark separation, with Zi — 1. In Eq. (jlip 
the gauge factors [zf, Zq] render the matrix element gauge- invariant; in the following 
the light-cone gauge A + — will be adopted, so that the gauge factors reduce to 
the identity. 

The matrix elements , and ultimately the DAs, are directly linked to the 

L z = component of the valence-quark wave function of the nucleon, by integrating 
out the transverse momenta of the constituent quarks. To the leading twist-three 
accuracy this amounts to writing the three-quark uud component of the proton 
state with positive helicity in the infinite momentum frame in terms of a single 
independent scalar function 3>(xi, X2, £3), where the variables x,, conjugate to the 
light-cone positions of the quark operators, are collinear momentum fractions of 
the proton longitudinal momentum carried by e ach quark, with < Xi < 1 and 
22i—i x i = 1 by momentum conservation. One has^ 



Wft,^ - -j&jl (n$=) * (i-i> ) *(*^>*3) 

WCki)u l (h)dH~k 3 )) - WCk^ihhHh)}] , (12) 




X 
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Table 1. The moments 0(*> m > n ) of the proton DA 
at different scales. The asymptotic values (AS) and 
the lattice results (LAT) at the scale Q 2 = 1 GeV 2 
are compared with the corresponding results from the 
present model calculation (Ref. 12), after evolution to 
Q 2 = 1 GeV 2 . 

(l,m,n) AS LAT Rcf. 12 



(0 0) 1 

(10 0) ± ~ 0.333 

(0 10) I ~ 0.333 

(0 1) I ~ 0.333 

(2 0) | ~ 0.143 

(0 2 0) i ~ 0.143 

(0 2) | ~ 0.143 

(110) J- ~ 0.095 

(10 1) J ~ 0.095 

(0 11) K- ~ 0.095 



1 


1 


0.3999(37)(139) 


0.340 


0.2986(11) (52) 


0.335 


0.3015(32)(106) 


0.326 


0.1816(64)(212) 


0.147 


0.1281(32)(106) 


0.144 


0.1311(113)(382) 


0.137 


0.1092(67)(219) 


0.098 


0.1091(41)(152) 


0.095 


0.0613(89)(319) 


0.093 



where /n is the proton decay constant. 

Correspondingly, with the LCWFs of the previous section one obtains 

/ 3 \ 



1~ 



N 



H 2( 27 r)3 



2(2tt) 3 
~7W 



i=l 



n ? 



2(2 7 r) 3 tf|J^k i xJ 

1/2] 

V({a:i,kii})H^(Ta,T), (13) 



where ip({xi, kij_}) is given by Eq. {§} and 



(T, |, T) = ^ J] -^== (-aia 2 a 3 + fcf *?a 3 - 2 ai fc 2 % L ) 



(14) 



with the following definitions: = (m + SjMo), N(ki) = [(to + .i^Afo) 2 + 
— ki x ~\~ ikiy : and = k{ X 




Fig. 4. The proton distribution amplitude. 
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In Fig. [4] the model results for the proton distribution amplitude $ are shown. 
The resulting shape is quite similar to that of the symmetric asymptotic DA <E> as = 
120x1X2^3. The results for the first and second moments cj)"' m > n > of $ with l+m+n < 
2 after evolution to Q 2 — 1 GeV 2 are shown in Table [T] Comparison with the 
analogous results from lattice estimated^ is quite nice. 



5. Transition distribution amplitudes 



The leading-twist decomposition of the general matrix element describing the tran- 
sition from a nucleon to a meson state, e.g. a 7r°, 

Tlp^ = {{^ {p-k)\ e* J ' fc Ua {z\n)u^(z2n)d^ {z^n) \p(p p , A))^ , (15) 

involves eight TDAsP^I The TDAs are dimensionless functions and depend on 
(xi,£, A 2 ), where the fraction of quark plus momentum Xi have support in 

il + A 2 M 2 



A 



-2£ 



- A 2 



(16) 



and the skewness variable £ describes the loss of plus momentum of the initial 
hadron in the proton-to-meson transition, i.e. 

A • n A+ >il „ 1 

= -, withP=- 

2P ■ n 2P+ ' 2 



(17) 





October 7, 2009 8:34 



proc'trento 



Probing the parton content of the nucleon 9 




Restricting ourselves to the case £ > 0, momentum conservation requires J2i x i = 
2£. The fields with positive momentum fractions, Xi > 0, describe creation of quarks, 
whereas those with negative fractions, Xi < 0, the absorption of antiquarks. This 
leads to define three distinct kinematical regions: the ERBL region for Xi > 0, and 
two DGLAP regions when x\ > 0, X2 > 0, X3 < 0, or x\ > 0, X2 < 0, x$ < 0. The 
names derive from the evolution equations which control the scale dependence of 
the TDAs in the different regions 

Results are presented in Figs. [5] and [6] for the TDAs in the ERBL region, cor- 
responding to probe the ip^q+qq and the contribution from the fluctuations of the 
nucleon in (jot ) and (A + 7r°) subsystems. The relative contribution of these com- 
ponents depends on the momentum transferred between the initial nucleon and the 
final pion as well as on the different spin configurations of the intermediate baryon. 
In particular, the A plays a special role in the case of the tensor TDAs Tf* which 
involve configurations with helicity ±3/2 (the only nonvanishing contribution to 
Tf"), while the interplay of the nucleon and A contributions with helicity ± 1/2 
determines the different shape of the vector (Vf*) and axial- vector (A? n ) TDAsM^ 



6. Concluding remarks 

In the light-cone description the nucleon state is decomposed in terms of A-parton 
Fock states with coefficients representing momentum light-cone wave functions of 
the N partons. Since the constituent quark models work so well phenomenologically, 
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in applications it is usually assumed that only the Fock components with a few 
partons have to be taken into account. One of such models has been studied in 
a series of papers to show that the parametrization of the LCWFs up to five- 
parton components is already sufficient to account for the electroweak form factors^ 



and spin densities^ of the nucleon, as well as the observed asymmetries due to 
transverse momentum dependence of parton distribut ions]^ and to give a useful 
insight into the quark generalized parton distributions jlfi U-7 | 18 | 2 | ^ s a f^-ther test 
of the model in this paper the nucleon distribution amplitudes and the nucleon-to- 
meson transition distribution amplitudes have been considered. 
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